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Abstract 

We discuss the scalar cosmological perturbations in a 3-brane world with a 5D bulk. 
We first show explicitly how the effective perturbed Einstein's equations on the brane 
(involving the Weyl fluid) are encoded into Mukohyama's master equation. We give the 
relation between Mukohyama's master variable and the perturbations of the Weyl fluid, 
we also discuss the relation between the former and the perturbations of matter and 
induced metric on the brane. We show that one can obtain a boundary condition on the 
brane for the master equation solely expressible in term of the master variable, in the 
case of a perfect fluid with adiabatic perturbations on a Randall-Sundrum (RS) or Dvali- 
Gabadadze-Porrati (DGP) brane. This provides an easy way to solve numerically for the 
evolution of the perturbations as well as should shed light on the various approximations 
done in the literature to deal with the Weyl degrees of freedom. 
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1 Introduction 



The phenomenology of brane world models has been the subject of intensive investigations 
in the last years. Its richness is in part due to the fact that those models can lead to 
modifications of gravity at small, but almost macroscopic, or even very large (cosmological) 
distances. The first case occurs in e.g. Arkani-Hamed-Dimopoulos-Dvali [Q] or Randall- 
Sundrum Q type of models, where the recovery of ordinary gravity at every-day scales is 
due to the geometry of the bulk space time, and gravity gets modified when Kaluza-Klein 
modes become excited. On the other hand, the Dvali-Gabadadze-Porrati (referred to as 
the DGP model in the following) brane worlds model leads to modification of gravity 
at very large (cosmological or larger) distances, in addition to possible small distance 
modification of different nature (see 1^-^) . In this model, the recovery of 4D gravity 
for usual distances does not rely on the geometry of the bulk, but is ensured by the 
effect of an induced Einstein-Hilbert term in the action, while the extra dimension can 
be non compact and have infinite volume 0. Because of those modifications of gravity, 
the cosmology of brane worlds can differ dramatically from standard cosmology and can 
potentially lead to various ways to test them. This is particularly true with the advent 
of precision cosmological measurement. Conversely the brane- world models can also lead 
to new scenarios for the primordial universe (like e.g. in |^) or its recent cosmological 
evolution (as in the DGP model, which has the ability to produce acceleration of the 
universe, as suggested by SNIa data without the need for a non zero cosmological 
constant [|lO|,|ll[ in a way currently in agreement with supernovae and Cosmic Microwave 
Background (CMB) data (see also [||,|l|])). The DGP model is also in itself an 

interesting playground to investigate the van Dam-Veltman-Zakharov discontinuity 

Measurement of the anisotropies of the Cosmic Microwave Background, and large scale 
galaxies or weak lensing surveys provide a unique way to test gravity at large scales (see 
e.g. ijl^) but also our comprehension of the physics of the primordial universe. As far as 
brane world models are concerned, this relies in particular on a better understanding of 
cosmological perturbations in these kind of models. The present article goes along this line 
and aim at reporting progresses in the study of brane world cosmological perturbations. 
The latter have already been investigated by many authors |p"7|-p^ with only limited 
results, as far as observable predictions are concerned (see e.g. pO|, p9[| ) . We will only 
consider here the case of scalar perturbations in a brane world model with a single 3-brane 
of codimension one. The bulk space-time will not be assumed to be stabilized, which is 
well suited for the DGP model or the Randall-Sundrum //model (presented in [||, that 
we will call the RS model in the following). In the latter case, it has been shown possible 
to obtain effective local perturbed Einstein's equations on the brane, which can be easily 
compared with usual perturbed 4D Einstein's equations Apart from the different time 
evolution of the background, and terms arising from the brane typical quadratic coupling 
between brane matter and induced metric [30|, the extra-dimension manifests itself by 

^This very fact can lead to new ways of addressing the cosmological constant problem 
^As far as the CMB is concerned, the model has only been tested in ^| assuming standard growth of 
cosmological perturbations, but non standard background evolution. One of the purpose of this work is to 
develop the formalism required to go beyond this approximation and test its validity. We will however, for 
the sake of simplicity, only give here the equations for a cosmological DGP brane which bulk space-time 
lies in its interior, which does not correspond to the case considered in jl^ (see po|). 
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source terms in the effective Einstein's equations which are coming from the bulk Weyl's 
tensor |31| and represent the imprint of bulk gravitational waves on the brane. Those 
equations do not close on the brane because of the lack of a local evolution equation 
for all the Weyl degrees of freedom^. It is however possible to close the system for large 
scale cosmological perturbations ||2^ , but this does not enable to compute the Sachs- Wolfe 
effect |2^. In order to obtain the evolution equations for cosmological perturbations on 
the brane, one should then solve the equations of motions for perturbations in the bulk. 
In this perspective, one may wish to use the work of Mukoyama who showed that the 
bulk scalar Einstein's equations are solved by a single master variable obeying a master 
equation p^. The master variable can be related to brane matter by junction conditions on 
the brane |2^. One of the purpose of the work reported here is to discuss the relationship 
between this approach and the previous one, and in addition the nature of the boundary 
condition for the master variable on the brane. This can ease the formulation of the 
problem if one wants to integrate numerically the equations of motion in the bulk, but 
also shed light on various approximations done in the literature to deal with the Weyl's 
degrees of freedom [^, ^ as well as on the well-poseness problem. Although most of 
this work can be applied to any kind of model involving a 3-brane of codimension one 
embedded in a bulk with a cosmological constant, we will, when considering realization of 
this setup, only deal with the case of RS and DGP model. 

The paper is organized as follows. In the rest of this section, we briefly introduce the 
models we are interested in and their geometry (subsections and |1.2D , and then review 
some properties of the background cosmological solutions for those models (subsection 1.3). 
We then introduce our notations for metric, matter and Weyl's perturbations (respectively 
section |2| and |3|), we also discuss causal properties of the master equation in relation with 
the issue of initial and boundary conditions (2.2), both for RS and DGP model. We then 
show explicitly, in section |^, how all the perturbed effective Einstein's equations on the 
brane, and conservation equations, are encoded into the master equation. In section ^, we 
obtain a boundary condition on the brane for the master variable for a perfect fluid on a 
RS or DGP brane, we also give the relation between the Weyl's fluid perturbations and 
the master variable. 



1.1 Models for brane worlds 

We consider a 3-brane embedded in a 5D space-time. The bulk space-time will be con- 
sidered empty with a, possibly non vanishing, cosmological constant A(5). The brane is 
defined as a hyper-surface X^{x^), where are bulk coordinates, and are coordi- 
nates along the brane world-volume^. The bulk metric induces through the brane 
^namely, for the Weyl's fluid anisotropic stress. 

^Throughout this article, we will adopt the following convention for indices: upper case Latin letters 
A, B, ... will denote 5D indices: 0, 1, 2, 3, 5; Greek letters /x, v, ... will denote 4D indices parallel to the brane: 
0, 1, 2, 3; lower case Latin letters from the middle of the alphabet: will denote space-like 3D indices 

parallel to the brane: 1, 2, 3. 
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embedding a metric on the brane, cahed the induced metric^ and given by 

ffit^ = d.X^d^X^gf^. (1) 

The gravitational part of the action of models under consideration can be taken as con- 
taining the 5D Einstein-Hilbert action given by 



S(5) 



-^/d^XV5(i)(%)-2A(5)), (2) 

'^(5) 



where A(5) is the bulk cosmological constant, i?(5) is the 5D Ricci scalar, and k^^^ is the 
inverse third power of the reduced 5D Planck mass. The models which we are interested 
in are (and have to be, in order to give an acceptable phenomenology) all sufficiently close 
to ordinary 4D GR in some region of their phase space. For that reason we also introduce 
a dimension-full quantity k'^^-^ which we set to be the inverse square of the 4D reduced 
Planck mass. 

To account for the brane, we add to the action (^) a term taking care of brane-localized 
fields given by 



5(4) = J d x^/g^^C, (3) 

where £ is a Lagrangian density that will be discussed below (we also implicitly include a 
suitable Gibbons Hawking term). 

It is always possible to choose a, so called Gaussian Normal (referred to as GN in the 
rest of this work) , coordinate system where the line element can be put in the form 

ds'^ = dy^ + gl^Jdx^'dx'', (4) 

and the brane is the hyper-surface defined by y = 0. In this coordinate system the 
gravitational equations of motion derived from the sum of action (^) and (^) are simply 
given by 

Gab = 4)S{y)S,,6^^6^s - A(5)5S, (5) 

(5) 

where G)^^ is the 5D Einstein's tensor, and S^i, is the effective energy momentum tensor 
for the brane. The latter is defined by the usual expression (where removing of total 
derivative terms is implicit) 




In the following we will refer to this tensor as the effective matter energy momentum 
tensor to distinguish it from real matter energy momentum tensor (see below). For most 
the calculations done in sections 0, H and 0, we will not need to specify the explicit form 



®We will in the following use an index to distinguish, when necessary, a tensor depending only on 
the induced metric (such as the induced metric itself g^^J , or the Ricci scalar computed with the induced 
metric and covariant derivatives compatible with the latter, -R(4)) from its bulk counterpart (such as the 
bulk metric or the bulk Ricci scalar) to which we put an index . 
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of £, which is model dependent. When deahng with existing models, we will consider two 
cases: 

(I) The Randall-Sundrum model // (RS model in the following) where the bulk is 
a slice of AdSs space-timeQ of radius r(5) given by 




and the brane is endowed with a tension (or cosmological constant) A(4) , and brane local- 
ized matter (that we will call real matter in the following) of Lagrangian density C(^]\jy In 
this case one has 



The brane tension and bulk cosmological constant fulfill the constraint 

A 4 (4) 

The effective brane energy momentum tensor is given by 

5,. = rW-A(4)#j, (6) 

where Tjuf^ is the real matter energy momentum tensor. The Kaluza-Klein (KK) spectrum 
of gravitons in the RS model contains a normalizable zero mode localized on the brane 
which dominates at large distance the gravitational exchange between sources. This gives 
rise to ordinary gravity at large distance, with a coupling given by 



K(5)A(4) 



(4)-—^ 

Below distances of order r(5), on the other hand gravity gets modified by massive KK 
excitations. 

(II) The brane-induced gravity model of Dvali-Gabadadze-Porrati (DGP model in the 
following) Q, where the bulk will be taken to be a slice of Minkowskis space-time^ (with 
vanishing cosmological constant and a vanishing Weyl's tensor) and an Einstein-Hilbert 
term computed with the induced metric on the brane is present in the brane action. The 
brane Lagrangian density is given by 

£ = £(M) - 

^'*(4) 



In this case Sny is given by 



^We will not discuss here the case of a brane in AdS-Schwarzschild, but a similar work is expected to 
apply to this case as well. 

*See for discussions of more complicated bulk geometry. 
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In the DGP model, the gravitational potential between static point-like sources behaves 
as in 4D (varying as the inverse distance between sources) for inter-source distances much 
smaller than a crossover distance Tc given by 

rc = (8) 

while it behaves as a 5D potential (varying as the inverse distance squared) for much larger 
distances. As mentioned in the previous subsection, one also expects to get small scales 
modifications of gravity, coming from quantum gravity effects, at distances smaller than 
^(5) @' This goes beyond the classical analysis of this paper and will not be discussed 
in this work. In the DGP model, there is no normalizable zero mode for the graviton 
excitation, and gravity, from a 4D point of view, is transported by massive KK modes. 
This can be understood as the source of the van Dam-Veltman-Zakharov discontinuity 
appearing in the DGP model in the linearized theory over a flat (Minkowski4) brane 
We will not discuss this issue here and refer the reader to Refs. [0,|3^-|35i. 



1.2 Some geometry of brane world 

Einstein's equations (^) lead to Israel's Junctions conditions which relate the jump of the 
extrinsic curvature tensor of the brane X^jy to whatever distributional source appears on 
the r.h.s. of (|5|) (here accounting for the brane localized fields or induced metric dependent 
terms). When a Gaussian normal coordinate system (^) is chosen, the extrinsic curvature 
tensor of a given y = constant hyper-surface is defined by 

K,u = \dygl^J. (9) 

Israel's Junction conditions derived from iM) simply read in this case 



-4) V - Isg^'j) > (10) 

where S = S^uq'^^^ is the trace of the effective energy momentum tensor of the brane, we 
have assumed a brane with Z2 symmetry, and the index (j,) means accordingly^ that the 
value of the component (here, as for any expression in the rest of this paper) K^y is taken 
on one side of the brane, namely in y = 0"*". The above equation (|To| ) can also be rewritten 

where K is defined hy K = K^^g'^^^ ■ 

Using a Gauss decomposition and the above Israel's conditions, one can then derive 
effective 4D Einstein's equations on the brane relating the brane intrinsic curvature to its 
effective energy momentum content |^^. In particular one gets the following equation |31] 

Gjt) = - Jff(t^A(5) + 4)"/.. - ^M-, (11) 



®When dealing with tensors, space-time indices will always be written without parenthesis, while other 
indices referring to properties of the quantity considered, such as , '''' , , will always be parenthe- 
sized. 
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where G^^} is the 4D Einstein tensor, li^y is given by 

TT _ ^ Q c"^ _L ^ Q Q _i_ ''^ c;" Q'^P (4) ^ q2 (4) 

and if^jy is defined as the hmiting value on the brane of the electric part of the bulk Weyl's 
tensor. In the Gaussian normal coordinate system (|^ it is simply given by 

£fj.u = fj.5u{b)j (12) 

where ^(j^ is the bulk Weyl's tensor. One sees that acts in the effective Einstein's 
equations ( [Tl[) as an external source with an energy momentum tensor that one can 
define as 

Tj^ = ——2—£/iu- (13) 

K(4) 

Following previous works, we will refer to this source as the Weyl's fluid on the brane. 
Because the trace of the Weyl's tensor vanishes, one sees easily from the definition (13) 
that in GN coordinate system (^) one has 

so that the Weyl's fluid shares some similarities with a radiation fluid. An other useful 
identity is the Codacci equation reading in a GN coordinate system 

where Z?^ is the covariant derivative compatible with the induced metric. Using the 
limiting value of the Einstein's equations (||) on the brane as well as the above equation 
one deduces that the effective brane energy momentum tensor on the brane is conserved 
with respect to the induced metric. From this and equations and (as well as from 
Bianchi identities in the latter case) one deduces that the brane matter energy momentum 
tensor is also conserved with respect to the induced metric in both models (I) and (II). 
One has then for the RS and DGP models 

£"^5^. = 0, (15) 



One also deduces from equation (11) and Bianchi identities that 



D'^S^, = ^jD'^n^., (17) 

which can be considered as a conservation equation for the Weyl's fluid. In general the 
energy momentum tensor of the latter is indeed not conserved with respect to the induced 
metric, since the right hand side of the above equation does not vanish, as will be seen 
below, we will however still refer to the equations deduced from (^) as conservation 
equations for the Weyl's fluid. 
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1.3 Cosmological background 



To deal with the background cosmological evolution, it is convenient to use a GN coordi- 
nate system (^) and to consider the following ansatz for the metric 

ds(5) = -n^{t, y)dt^ + a2(t, y)6ijdx'dx^ + dy"^, (18) 

where the 3D metric 6ij is a flat euclidean metric (we will only consider in this paper the 
case of a spatially flat universe). One can further choose a time parametrization such that 
the function n is set to 1 on the brane (n(ft) = 1), in which case the induced metric ([l|) is 
simply given by 

'^5^4) = —dt'^ + a^fj^-^5ijdx^dx^ 

(where a(6) = a(t, y = 0)) and is of FLRW form. Considering comoving observers to sit at 
fixed comoving coordinates on the brane, t is then simply the cosmological time on the 
brane. Accordingly with the symmetries of ( [T^ ) the effective energy momentum tensor of 
the brane is taken of the form 



S^ = 5{y)<Xmg{-p,P,P,P). 

We point out that in the above expression p and P are not to be understood as real matter 
energy density and pressure, that will be denoted /0(m) aiid P(m)-, but as their effective 
matter counterparts. Their relation to real matter energy density and pressure is model 
dependent. For example in the case of RS model one has 

P = P{M) + A(4) • 

The expressions for p and P as functions of P(M) ^-iid P[m) in the RS and DGP models 
are given in appendix |A.2| . The junction conditions (|lO|) are given for the background by 

P (19) 



= ^(3P + 2p) (20) 
n{b) 6 

where a prime denotes a derivation with respect to y. The behavior of the scale factor on 
the brane can be obtained solving the equations ( pT|) and (|l^) which for the background 
reduce to the two independent equations valid for a bulk with a vanishing background 
Weyl's tensor (we will always assume this to be the case in the rest of this paper) 

H' = + (21) 

p = -W{P + p), (22) 

where a dot means a derivation with respect to time t, and H is the Hubble parameter on 
the brane given hy H = a(fe)/a(fe). 
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One finds in the Randall- Sundrum model, using the above equations, a transition 



1 36, 37] from an early non standard phase to a late time standard cosmological phase 
which dynamics is given at leading order by usual Friedmann's equations 

2 

■J 



t2 _ ^^(4) 



-P{M), 

P{M) = -^H{P(^M) + P(M))- (23) 

The transition happens approximately when the Hubble radius becomes larger than 
the threshold r(5) . 

For the DGP model, one has a transition from an early standard cosmological phase. 



governed by equations (p3|), to a late time non standard phase |1C|. The transition happens 
approximately when crosses the length scale Vc- A conservative bound on Vc (in 
agreement with other bounds, see ^) is then given by demanding that rc should be 
larger than (or of the order of) today's Hubble radius. As a consequence the parameter 
T defined by 

is always smaller than (or of the order of) one in the past. 

An explicit expression of the bulk metric can also be obtained in the coordinate system 
( p^ p7[ | , we quote in appendix ^ some useful identities concerning the solutions for the 
background metric both in the bulk and on the brane. 



2 Scalar metric perturbations in the bulk 



In this section, we first (subsection |2.1| ) introduce our notations (mostly the same as 
in |2^) for scalar metric perturbations and gauge invariant variables. We then remind, in 
subsection 2.2, the expression of Mukohyama's master equation |24| as well as the relation 
between the master variable and gauge invariant variables. Eventually we discuss there 
some properties of the Master equation in relation with issues related to initial conditions 
in the bulk and boundary conditions on the brane for the DGP and RS models. 



2.1 Metric perturbations and gauge invariant variables 

The five dimensional perturbed metric element is taken to be 

/ -n2(l + 2A) a'^B\i nA 



gAB 



V 



(1 + 2n)6ij + 2E\,j 




(25) 



where we have kept only scalar perturbations, and |j denotes a differentiation with respect 
to the comoving coordinate x*. A scalar coordinate transformation is defined as 



t 



t + 6t, 



(26) 
(27) 
(28) 
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and induces the following transformation on the scalar quantities 



A 

n 

B 

By 

E 



A- 

n 



6t 

a 
a 



n 



6t 



n 



n 



n 



6t 



B + -^St 

By — 6x' - 
E-6x, 



Ay + n6t' 5y, 



—Sy, 
a 

- 6x, 

X^y, 



n 



A, 



yy 



A, 



yy 



6y'. 



In the following a gauge transformation involving a zero Sy will be denoted as a 4^- 
gauge transformation as opposed to a gauge transformation involving a non zero 6y which 
will be referred to as a 5D-gauge transformation. For further reference we note that the 
(perturbed) brane position can be defined locally hy y = We will also use the 3D 

spatially gauge invariant variables |20| 

a = 
ay = 

transforming as 



-B + E, 

-By + E' , 



a 



a 



-^St 
1 



cr,. 



under the coordinate transformation (12' 
under 4D-gauge transformations by 



We also define the variable ^ and $ invariant 




aa 



tCT. 



From equations (26-28), one sees that the most general (scalar) 5D gauge transformation 
involves 3 unknown functions, implying that one can define out of the 7 perturbation 
variables appearing in (|2^), four independent variables invariant under a 5D-gauge trans- 



formation |^,|23]. We take the latter to be given by A, Ay, Ayy,R, expressed as 

A = 



A. 



yy 

n 




n 



a, 



Ayy + (a^cr. ^' 
Tc + aa a-u kCt. 



10 



2.2 Master equation 



As shown by Mukohyama |24|, the hnearized scalar Einstein's equations in a maximally 
symmetric bulk can be conveniently solved introducing a master variable which obeys 
a PDE in the bulk, the master equation. The latter, when Q has a non trivial dependence 
in the comoving coordinates rr*, reads in a GN coordinate system (18) 

A\ nn fnn'V , , 




where A is defined by A = b^^djdi. In the rest of this article, we will implicitly consider 
all the perturbations as Fourier transformed with respect to the x*s, in order to do a mode 



by mode analysis. In particular (29) can be rewritten as 

VtSi = 0, 

where Pa is a second order hyperbolic differential operator acting on y and t dependent 
functions (in the GN system), and A is understood to be replaced by —k'^, where k is 



the comoving momentum. Equation ( |29[) is then only valid when k does not vanish |24|, 



which is the only case of interest as far as cosmological perturbations are concerned. The 



gauge invariant metric perturbations are then related to by 1 2C , 2^ 

— (in" + irn + -"\. 

6a V 6 n-^ n 

Ay = — fn' - -n] , (31) 

an \ n J 



A = -±(2n" + A.n+'-^n-^n--n'], (30) 



Ayy = ±ln" + 4^n-'-^n-2^n-2-n') , (32) 



1 / „ 2 ■■ Af- " 

n = ^(n"-l,n + ^n + ^n + -n'). (33) 

ba \ 3 ri'^ n / 

The master equation enables a priory to solve for the evolution of brane world cosmological 
perturbations, once initial conditions in the bulk and a suitable boundary condition on 
the brane are provided, as we now discuss. 

Let us first deal with the DGP model, where the bulk is simply a slice of 5D Minkowski. 
In this case, the metric ( |T8|) can be put in the form 

ds^ = -dXdT + X^dx'dx^ij (34) 

where X is given by a(t, y) and T is a function of the Gaussian normal coordinates y and t 



which can be easily obtained e.g. by the coordinate change given in |39|. The background 
trajectory of the brane X(^)(t) and T'(b)(t) is then given by 

T(b) = a^^J, 

with a(fc) obtained by solving equations (^) and (p2|). In the coordinate system (X, T), 
equation (|9|) reads 

= -dxdrn + + ^ = 0. (35) 
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where the hyperbohc nature of 2?a is manifest. The characteristics of Pa are simply the 
null lines X = constant, T = constant. If one assumes that d(;,) is always positive, then 
is an increasing function of time, in which case the trajectory of the brane can be 
schematized as shown on figure ||. The brane background trajectory is then non char- 
acteristic. Let us now assume that a given initial time is chosen along the brane, from 
which one wishes to evolve cosmological perturbations. One then needs to specify initial 
data in the bulk , e.g. as shown in figure ||, as well as a proper boundary condition along 
the brane. We will not discuss here how the bulk initial conditions are provided^'^; we 
rather aim, as far as the issue of initial and boundary data are concerned, at clarifying 
the expression of the boundary condition on the brane. The latter can be imposed along 
the background brane trajectory, and will be given and discussed more explicitly in section 
5.2'^. We further note here that the bulk initial data may involve in general character- 



istics initial data0 (meaning here data specified on a characteristic X = constant line). 
Eventually, we stress that the coordinate system (X, T, x*) is singular in X = 0. One can 
further transform the metric (^4|) to a flat Minkowski metric of the form 

ds^ = -UY+dY- + {dY^f + {dY^f + {dY^f, (36) 

where Y~^ and Y~ are light cone coordinates and one has X = Y^ . In this coordinate 
system, the master equation reads 

Pa = r^^^dAdB^ + ^ = 0, 



where tiab is the metric appearing in (|3q). Using this coordinate system, one can easily see 
that a brane with a Big Bang is a conoidal hyper-surface with a vertex that can be taken 
to be the origin. The brane is then inside the light cone of the origin which represents the 
Big Bang singularity. The latter also extends along the null line y~ = on the light cone 
of the origin, to which the brane is tangent in strong analogy to what is happening when 
the bulk is a slice of AdSs One notes that the master equation is singular along this 
line, where the scalar-vector-tensor splitting of Minkowskis gravitational perturbations 
becomes ill-definedP| This is not without importance if one wants to set initial data close 
to the initial singularity. 

In the case of the RS model, the bulk is a slice of AdS^ space-time, and the metric 
( [Tsl) can be put in the form 

ds^ = ^ (^-dr^ + dz"^ + Tjijdx'dx^'^ 



where one has z = r(^^-j/a{t, y) and r is a function of the GN coordinates y and t that can 



be found using the work of |^l[. The brane trajectory is given by |42| 



a 



^'^The latter can be thought of as being given by some early universe physics. 

^^We also discuss more in appendix the issue of bulk initial data in relation with the boundary 
condition found in this section. 

^■^as opposed to Cauchy type of data specified on non characteristic curves. Such characteristics initial 



data can also be considered as a model dependent boundary condition in the bulk 

^■^This is because this splitting is defined by a particular choice of coordinates on Minkowskis, e.g defined 
by equation (^^, which breaks down along the = X = line 



12 




Figure 1: Schematic representation of the bulk space-time with the brane trajectory in 
the characteristic X,T coordinates (of equation (p^)) for the DGP model. The left figure 
is the case with a Big Bang (with e.g. a radiation dominated universe), while the right 
one corresponds to a case with inflation. The X = Q line is a coordinate singularity. The 
gray (cyan) part is cutoff and the complete bulk space-time is made out of the remaining 
part, glued to a copy of itself along the brane. If a given initial instant is chosen along the 
brane from which one wishes to evolve cosmological perturbations, one needs to specify a 
boundary condition along the brane (AB), and initial data in the bulk (e.g along (AD)) 
that may be in part characteristic data (as along (AEF)). The future of the brane events 
with cosmic time larger than initial time lies above the dotted line (AC). 



1 

a{b) 



\ 



1 + ^2 !k 



In the coordinate system {z,t), equation (^) reads |^8| 

= din - a^f) + ^ + ^ + Af) = 0. 

Using the null coordinates X = t + z and T = t — z, this can be further rewritten 

417 



Va = —Oxori'l H — ^ h TTT 

A X 1 ^ ^ _ri. ^ (X - r)2 



+ AH = 0. 



(37) 



We also see clearly the hyperbolic nature of The characteristics curves of the latter 
are the null lines X = constant and T = constant, and the brane trajectory is also 
non characteristic. We will not discuss here in more details this case since it is strongly 
analogous to the previous one. 
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3 Matter and Weyl perturbations on the brane 



In order to obtain a boundary condition on the brane and to recover perturbed effective 
Einstein's equations from the master equation, one needs relations between matter, Weyl's 
fluid, and the master variable. These can be obtained by using limiting values on the 
brane of derivatives of bulk metric perturbations which are related to effective matter and 
Weyl's fluid perturbations (that we define in the next subsection) through Israel's junction 
conditions (|To|) and the Weyl's fluid definitions (|l2[]T3|). The corresponding expression, in 



a convenient gauge defined in subsection |3.2| , are reminded in subsection |3 



3.1 Perturbations of effective matter, real matter and Weyl's fluid energy- 
momentum tensor on the brane 

Following standard cosmological perturbation theory, we can always decompose the scalar 
part of the perturbations of the brane effective energy momentum tensor as 

5S^ = -6p, (38) 
6S^ = <5(Z|„ (39) 

1 

— ( 
3 



<55j = 6P5} + a;. - -6) A ] 6Tr, (40) 



where A*- is defined by 



so that one has A = A^. Similarly the scalar real matter perturbations are defined from 
perturbations of the real matter energy momentum tensor and read 

'^^0%/) = -^P{M), (41) 
6Tj^^M) = ^mm, (42) 

6T;^m) = 5P(M)5j+(A;.-i<5;.A)57r(M). (43) 

The relation between the effective matter perturbations (38-^) and the real matter per- 
turbations (^^43) in the two different models 1/ and II/ mentioned in the introduction 
can easily be read off from equations (|6|) and (^) respectively. They will be given in section 
5.2 . Similarly, the Weyl's fluid scalar perturbations can be decomposed as |21, 2^] 



5£lj = kI^)6p(s), (44) 
^£•0 = -fq^^^Sq(^^^\j, (45) 

58] = -4)(5J^,)5j + (Aj-i5^A)5vr(,)). (46) 

Prom equation (|l^ ) one can deduce that the Weyl's fluid perturbation obeys the 'equation 
of state' of radiation, namely one has 

^P{S) = \5p{e)- (47) 
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3.2 Gauge choices 



From equation (l2g), one sees easily that one can simultaneously choose a gauge where 
are vanishing and the brane sits at y = |2^. This gauge choice put the 
perturbed metric in the GN form (^). This however does not completely fix the gauge, 
but one can further make y — independent shifts of x* and t without spoiling the gauge 
conditions enumerated above; this gauge freedom is simply the usual 4D gauge freedom, 
and one can use is to go to the usual (4D) longitudinal gauge on the brane. We will refer to 
such a gauge choice as a Gaussian Normal Longitudinal gauge (denoted as a GNL gauge) 
in the following. To summarize, the GNL gauge is defined by the conditions 

Ay = 0, (48) 
By = 0, (49) 

Ayy = 0, (50) 

= 0, (51) 
cT(fe) = 0, (52) 

where the last equality defines the 4D longitudinal gauge on the brane, and we have used 
bars to differentiate quantities in the GNL gauge from quantities in arbitrary gauge (e.g. 
^ vs. A). In the rest of this paper, we will always work in the GNL gauge when dealing 
with quantities of which we take limiting values on the brane, and in order to alleviate the 
notations we will drop the bars. In the GNL gauge, one has the following identities valid 
in the bulk (say in the vicinity of the brane) 

Ty = E', (53) 
a = -B + E, (54) 




A = A--{ + —a'E', (55) 

n \ n I 

Ay = ^ '- + ^ '--2^a'a, 56 

Ayy = {a^Ej, (57) 

ad 

-^a, (58) 



and on the brane 



7^ = 


n + aa'E' - 


{A = 

in = 






^ + !!^a''E') 
n ) 


Ay = 


{a^E'y a" 
+ — 

n n 



/ (fe) 

a^E" + 2a'aE') , 
-^ + aa'ii;')(fe). 
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We underline that the gauge condition ( |52[ ) does not imply the vanishing of y— derivatives 
of a on the brane, that have been kept above. 



3.3 Junction conditions for effective matter and Weyl's fluid 

The Israel junction conditions (|^) in the GNL gauge lead to the following expressions for 



the 7/— derivatives of the metric perturbations on the brane |20| 

,2 



= ^i36P + 26p), (59) 
7^'(,) = ^4)(A57r-M, (60) 

2 



= -If^y^- (62) 



One has also 



"(fe) 

Similarly, one can get expressions relating the second y— derivatives on the brane of metric 
perturbations to Weyl's fluids (and effective matter) perturbations. They read in the 
GNL gauge [20| (those expressions are obtained using the expression of the Weyl's tensor 



perturbations (4^^6|), the definition ( [T2|) and some bulk perturbed Einstein's equations) 

"(4) ('5%) + ^*P(£)) = + 

- -5|{-"-(4-^)-'}„. <-) 

= \e" + 2-E'] , (65) 
^ ' I « J (6) 

nf^-^Sp^S^ = l^3(^R" + 2^R'^ + AE" + 2^AE'^ . (66) 

One can consider the above expression ( |59[]6^ ) and (|6^|66[) as definitions of effective mat- 
ter energy momentum perturbations (e.g. 5q) or Weyl's fluid perturbations (e.g. Sq(^£)) 
respectively in terms of first (e.g. B') and second derivatives (e.g. B") transverse to the 
brane of bulk metric perturbations. There is indeed a one to one correspondence between 
the former and the latter, and we will in the following call all equations (pS|-|66|) junction 
conditions (respectively for effective matter or Weyl's fluid perturbations). Following this 
line, we have not simplified equation (^) using the equation of state (^7|), to keep the 
correspondence manifest. We will do so until the end of next section. 
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4 Linearized effective Einstein's and conservation equations 
on the brane from the master equation 



With obvious notations, the perturbed 4D effective Einstein's equations (11) reads 



(67) 



where one can insert the definitions (38-^) and (44-46) given in the previous section to 
get expressions relating the perturbed induced metric on the brane to effective matter 
and Weyl's fiuid perturbations. Those expressions can then be easily compared to usual 
4D cosmological perturbation equations p^. The latter are given in appendix |B.1|. The 



components of the linearized effective Einstein's equations are gathered in appendix |B.2 
together with the ones of linearized conservation equations ( [l5| ) and (^) . 

We would like here to show how those last equations can be recovered from the Master 
equation ( [29| ) and the junction conditions for effective matter and Weyl's fluid given in 
previous subsection. The outline of our derivation is to built 'constraints' valid everywhere 
in the bulk on the gauge invariant variables A, Ay, Ayy and R (as well as on their y— and t— 
derivatives) out of the master equations and the definitions ([30^331) . The latter constraints 
can be thought of as reconstruction of linearized 5D Einstein's equations in the bulk out 
of the master equation, in confirmation of Mukohyama's work |^^. We then take limiting 
values on the brane of those constraints, where we replace the gauge invariant variables by 
their expressions as functions of effective matter, Weyl's fluid and induced metric that can 
be obtained through the junction conditions. Those will yield the sought for equations. 
Let us now proceed. 

First, using equations (55-|5^), the derivatives with respect to y of those equations, the 
junction conditions for perturbations (given in subsection |3.3| ) and background (equations 
([l9|-^)), as well as the equation of motion for the background (gathered in appendix |A. 21) , 
we get the following expressions relating limiting values on the brane of gauge invariant 
variables (and their y— derivatives) to effective matter, Weyl's fluid and induced metric 
perturbations 



A 



A. 



Vib) 



A 



yy{b) 



7^, 



{b) 



(b) 



$ - -^a" (3P + 2p) 67r 



{-'t(5)'5g-'«(5)«M^)'}(^) , 



-^(6) + ^«(fe)P'57r, 



(68) 

(69) 
(70) 

(71) 



\5) 



(5). 



{36P + 26 p) + K^5)54 + -^a (a<57r) 



6 



(72) 



J J (b) 



A' 



y{b) 



(5p + 6P) [aiadny +6q] 
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7^'(,) = |-^5p + -|)A5^ + ^d(a<57r)- + ^4)<5(? 

''(')-4)pa25vr(,)| . (74) 



6 



The value on the brane of A'yy can however not be computed at this stage as can be 
seen from equation (|57|). It would indeed require to know E"l^^ which is not given by the 
junction conditions. Let us now derive three constraints on the gauge invariant variables 
A, Ay, Ayy, R aud their first derivatives A',Ay,R', which limiting value on the brane we 
know from the above expressions. 

Adding equations (^), (|3^) and (^3|) one first gets that everywhere in the bulk p0|,p7| 



A + Ayy + 7^ = 0, (75) 

which is the first constraint we will use latter. 

We now consider the five equations obtained taking first derivatives of equations (|30[), 
(|3l|), (|3^) with respect to y and t. They read symbolically 





h 


n- 




6 






n 




n 


n- 




6 


rfi 




n 



A' = -{^ (2^" + \^+ ^-^n -^n--n')[ , (77) 



6a V n 

Ay = , (78) 



an \ n 

These five equations, together with equations (|30|), (|3l|), (^), the master equation ( p9[ ) and 
its derivative with respect to y and t can be considered as forming a system {Sq) of eleven 
equations for the ten unknown functions Z// = {Q., VL' , $7", $7', il", ', il "', Q'"}. This 
means that one should be able to obtain a new constraint from it, valid everywhere in the 
bulk. Namely one gets 

Ayy + -Ayy + 2A+ 2- A + Tl'^^ + —Ay + - A' = - — + — 

a a " 2a ^ 2 y 2a V 6 n 

where the right hand side of the above equation vanishes, when taken into account the 
equations of motion for the background leading to ( |124 ). If we subtract to this equation 
two times the time derivative of the constraint (l75|) , we obtain the equation 

- Ayy + "^Ayy + 2^i + u'Ay + Ay + ^^A'y - 211 = Q , (8l) 
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which will be used latter. Removing equation ( pO| ) from the system (Sq), we get a system 
(S) of ten equations for the ten unknowns U with second members given by gauge invariant 
variables and their derivatives. This system has a non vanishing determinant in the bulk 
and can be inverted to give Q and its derivatives as a function of the gauge invariant 
variables. E.g. one obtains the following expression for $7 (which has been simplified using 
the equations of motion of the background metric) 

A^Jl = A(5) a^a ^— 2d(ft) + a'Ayy — aA' 



3a3 



+aa' (4A(i + Ayy) + 6dJ;,)(i + 2Ayy) - 3h^b)a'Ay 
+6ad(^) (2d(b)Ay + a A' + 2a A 



+3a(i(fc)a (id^^) - ^yy) + '^aA'y + a'Ay^ , (82) 

which shows in particular that J7 is 5D-gauge invariant. 

In a similar manner to the derivation of equation (pW , one deduces from the equation 



the following relation 

na' — an' ~ ~, 2a'n + an' ~ a ~ 1 ~ ~, $7' / A5 n" 

A -A' + Ayy Ay A -271' = — + — 

an an 2an 2n ^ 2o V 6 n 

which leads to the constraint 

naJ-a^^ _ ^ ^ 2a!n^ _ a^ _ 1 _ ^ ^ 
an an 2an ^ 2n ^ 



Taking the limiting values of the constraints (|7^), ( |8l| ) and ( |83[) on the brane, and 
replacing in it the gauge independent variable by their expression as a function of effec- 
tive matter and Weyl's fluid given by equations ( |68[]7^ ) one obtains the following three 
equations 

U-^ + a\l)67r^£)-'^a^{W + p)67r] = (84) 
I J (fe) 

i^24> + 2^<^+'^6qp + Kf^^6q^e)^ = (85) 

lsq + 3-6q + 5P + {P + p)<^ + -A6ir\ = 0. (86) 
la 3 J (b) 

Those three equations are respectively the trace free part of the ij component of the 
perturbed effective Einstein equation (|67|), the Oi component of the same equation, and 
the conservation equation for the effective matter momentum 6q (that is to say the i 
component of the perturbed equation (|l5|)). Those equations are derived in appendix BJ. 



from equations (^) and (^). We mention here that equations ( |8^ ) and (p5|) can be used 
to obtain an expression of (57r(^) and (5g(£:) as a function of effective matter perturbations. 
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The latter can then be used in equations (38-74) to express the left-hand side of those 



equations solely in terms of the effective matter and induced metric perturbations. 

We would like now to derive the remaining perturbed effective Einstein's and conserva- 
tion equations. Those imply constraints involving second derivatives of the gauge invariant 
variables. To be more precise, we will need the limiting values on the brane of A' , Ay, Ay 
as well as R",A". The first three limiting values are simply given by time derivatives of 
equations (|^), ( |69| ) and (^). However, computing the limiting values of R",A" on the 
brane require knowing E'" on the brane, as can be seen from the expressions (^5|) and 
(|5^. This can be obtained using the first constraint (75), leading to 



{a'E'" + ia'aE" + {2aayE'] 



where the last equality comes from (57). Solving for E"!^-^ one obtains (using again the 



junctions conditions of section 3.3 as well as the background equation of motion of appendix 



A. 21) 



6tt — adK^g^Jvr 



^(5) 

72 

(5) ,.2 2 



|24a2A(5) + Kj5)a2 (llp"^ + 12Pp + OP^) + 12A} ^vr 



+^^i4)^"{p-P)SH£)\ , (87) 

^ J (6) 

which allows to compute R"j^^ and ^"b)- The next step of our derivation is to extract the 
fourth derivatives Q " i} "", O '", Q"" from the expression of Ay and Ay as a functions 

of Q and its derivatives (that one gets from equation (^)) , and from the second derivatives 
dtdt, dydt, dydy of the master equation (p9[). Having done so one can insert those fourth 
derivatives in the expressions of A' , R" and A" as a function of and its derivatives. This 
lead to an expression for each of the quantities A', R" and A" which does only contain 
derivatives of of order lower or equal to three. One can then use, as in the derivation of 
constraints (|75|), (|8l| ) and (^) the system {S) to obtain an expression in the bulk relating 
each of the functions A\ R" and A" to Ay, Ay, A ,A', Ay, A'y,Ayy, A, Ay, Ayy as well as to 
the background metric. The last step is to take the limiting values of those expressions 
on the brane, and to replace in it all the gauge independent variables, as well as their 
derivatives, by their expressions as a function of brane effective matter and Weyl's fluid. 
The expression for A' leads to 

{5p + ^-{5p + 5P) + ^6q-2,{P + p)i!\ =0, 
La } (b) 

which is the component of the conservation equation (^), or in other words the effective 
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matter energy density conservation equation. The expression for R" leads to 

which is the 00 component of the perturbed effective Einstein's equation (|67D. While the 
computing on the brane A" + 2R" leads to 



-6-^ + -^A($-vI/)| , (90) 



which is the trace of the ij components of effective Einstein's equation (|6^). The last 
two equations are the equations coming from perturbing the Weyl's fluid conservation 



equation They can then be derived from equations (|4|), (|5l), (H), (H), (H), (|0D, 
(as a consequence of Bianchi identities) and read 

{^P{e)+\{^P{s) + ^Pie))+^2^m] ^ = 0' (91) 



and 

A 



Uq + 2,-5q^e) + SPis) + ^A57r(£) 1 =\^{P + p) U-5q + Adn -6p)\ . (92) 
L a -3 J (6) \ a y J 

We have thus showed at that point that all the scalar components of the perturbed effective 
Einstein's equations ( [Tl| ) and conservation equations (|l5|) and (^) can be derived from the 
Master equations and the junction conditions. We now discuss the form of the boundary 
condition for Q on the brane. 



5 Boundary condition on the brane from real matter equa- 
tion of state 

We first note that inserting the expressions (|68|-[7^) (as well as using equation (^6|)) in the 
limiting value on the brane of equation (|8^), one gets the following expression for A^0(5) 

A20(b) = |_i8a^d(^'f + ^$^ +6a3A^'--^oVp| ■ (93) 
Using similarly the expression for ^1' obtained from (S) we get 



An[^-^ = l-^pAn + nfr^^a^6p-3Kfr^^a'^a5q-Kf^^a^A6TT\ . (94) 



) (h) 



Similar equations were obtained by Mukohyama [28| and used by him get an integro- 
differential equation for the effective matter and induced metric perturbation on the brane 



21 



One may be tempted to consider them as a way to fix a boundary condition for the 
differential operator Pa- This is however not the case since the right hand side of those 
equations is not known as a function of time. Moreover, these equations do not allow to 
easily disentangle the extra physical conditions that have to be imposed to get a boundary 
condition, nor to discuss the well-poseness issue. To obtain the latter we wish to express 
the effective matter and induced metric perturbations as functions of the master variable, 
as we now do, discussing first the case of Weyl's fluid perturbations. 



5.1 Relation between Weyl's fluid, eff'ective matter perturbations and 
the master variable 

5.1.1 Weyl's fluid 

From the equations ( |3^ and ( |70|) one easily finds an expression for the Weyl's fluid 
anisotropic stress 5'/r(£:) given by 

'^^(^) = 7rA^|3^^-3^i^-^J^-^4)(P + p)f^'l . (95) 
o'^(4)"(f,) L a u z, J 

The energy density perturbation of the Weyl's fluid can similarly be derived from equations 
@) and (|8|), one gets 

The perturbation of the Weyl's fluid pressure is then given by equation (^7[), whereas 
can be obtained using the conservation equation (^) and reads 

Those expressions can be used to compute the perturbations of the Weyl's fluid energy 
momentum tensor once il. is known as function of time, and e.g. compute the Sachs- Wolfe 
effect ||2^|3. They can also be compared with similar expressions that will be obtained 
below for effective (or real) matter and induced metric perturbations in order to test 
various approximations done in the literature to deal with the Weyl's fluid (like in [p6|). 



5.1.2 Effective matter 

In analogy to the above expression for the Weyl's fluid perturbations, we wish to obtain 
expressions for the effective matter energy momentum perturbations and induced metric 
on the brane as functions of the master variable J7 and its derivatives. To achieve this, one 
can think to use the system (5), which can also be considered as a system (5') in the 10 
unknown functions 5p,5P,5q,6Tr,6TT,6Tf,^,^,^,^ once one replaces the gauge invariant 
variables (and their relevant derivatives) by their limiting values on the brane as functions 

^''The Sachs- Wolfe effect can als o be comp uted di rect ly, under the same hypothesis, using the expressions 
for $((,) and given in ( |l07| ), ( |ll3| ) and (|lo|), (|ll^). 
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of the effective matter and induced metric perturbation. This system (S') is however not 
invertible, but one can obtain from it the fohowing expressions 

,^2 




A 



(5) 



n + 6-n 

a 



3n + ^ (3P + 4p) n' 



a-" (3P + 2p) 6tt 



(b) 



6p 







^(5) 






2 




An 


a 

+ 3- 

a 


6a(f,) 



n + 3-n + ^pn' + 

a 2 



a pdir 



(98) 



(99) 



(6) 



an' 

-IS— — +6aA(57r 



18da (aJvr)' 



(fe) 



6a 



{b) 



{3P + 2p)n- -^n ' - 6a^ {aSiry 



(100) 



(101) 



6P 



1 



6o(b) 

+pn 



{P + p) 



2A A 



(5) 



n 



{b) 

a 



3P + - (4p + 6P)\n 



\5) 



{P + p) (3P + 4p) 



a 



^^n' + -^n-' 

'^(S) '^(5) 



4— + 4A(5) + 



"(5) 



(9P^ + 9Pp + ip' 



6TT + 6{a^6TTy 



(102) 



(b) 



One can verify by inserting those expressions in equations ( |8^ ) , (H) , ( |86| ) , ( |88[ ) , (|8^) , ( |90| 



(|9lD and (|92|) that those equations, namely the perturbed effective Einstein's equations 
and conservation equations, are identically satisfied whatever n. This means that the 
latter do not contain more informations than the above equations. The fact that (S') is 
not invertible could have been expected counting the number of extra equations one needs 
to get a closed system (or a boundary condition) as we now explain. 

Let us first recall what is happening for usual 4D cosmological perturbations (the 
linearized 4D Einstein's equations and conservation equations are reminded in appendix 
B.l ). In this case if one considers only one matter fluid as a source in the Einstein's 



equations, and using a decomposition of its energy momentum tensor of the form given in 
equations (^43), one gets a system of equations from the perturbed Einstein's equations 
which does not close. To close the system and solve for the evolution of cosmological 
perturbations once initial conditions are provided, one needs to add more equations, e.g. 
equations of state for the matter fluid. The simplest case considered in general is the case 
of a perfect fluid, which implies 



with adiabatic perturbations obeying 



(103) 



(104) 
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where is the sound velocity defined by 



P(M) 



(105) 



The two conditions (103) and ( |104| ) are sufficient to close the system of perturbed Einstein 
equations, leading to the evolution equation for the gravitational potential $ (defined as 
in equation ( 130[ )) 



^ + (4 + 3c2)^6 + 



d2 



2- + — (l + 3cf 

a a 



(106) 



and all the other relevant quantities can be deduced from the knowledge of $ as a function 
of time. 

Turning back to the brane perturbations, if one would be able to get an expression for 
^"^{M) ) ^P{M) and 5P(^M) as a function of 17 and its derivatives on the brane (that is to say if 
(5') would have been invertible), then imposing the two conditions (|103| ) and (|104| ) would 
lead to two independent boundary conditions on the brane, which is too much, once initial 
conditions are supplied on a suitable hyper-surface in the bulk. From the above discussion 
one can however expect as well to have to use two extra equations (e.g. equations of state 
on the matter fluid on the brane) in order to get a solvable problem for the evolution of 
cosmological perturbations on the brane. The first of those extra equations can be though 
of as enabling to express the matter fluid energy momentum perturbations as a function 
of the master variable and its derivative (namely to invert (5')), while the other leads to 
a single boundary condition on the brane. This will be done explicitly below for the two 



models described in section 1.1 



5.2 Boundary condition for real perfect fluid with adiabatic perturba- 
tions 



In the case of the RS model (model (I) in section |1 . 1| ) , one has simply the following 
relations between effective and real matter energy momentum perturbations 

= ^P(M) , 
6P = 6P^M), 
dq = Sqi^M), 
6iT = <57r(M) , 

as can be seen from equation (^). We now specialize to the case where the real matter 
is a perfect fluid. In this case the anisotropic stress (^'/r(jv/) vanishes and one gets from 
equations (p8-102) 
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2 

+pn-'^{p + p){2,p + Ap)n' + + -^0.-') . (Ill) 

^ '«(5)« '^(5) J(fe) 

These expressions allow to get a boundary condition for the operator Da on the brane, 
using an extra physical condition on the real matter and induced metric perturbations. 
In the case of adiabatic perturbations, the latter is given by equation ( [1041 ). If we insert 
in ( |104|) the expressions for 5p(Af) and 5P(^m) given by equations ( 1091) and ( |lllD we get a 
boundary condition for Da which is expressible by 

i^fl5(S^)(fe)+G,i5(f^')(fe) =0, (112) 
where F^is and Grs are polynomials of the cosmic time-derivative, dt, with cosmic time 



dependent coefficients which are known from equations ( |105|) , ( [109D and (111) as well as 
from the background solution. 

We now turn to the DGP model. In this case, one gets from equation (0) the following 
expressions for the effective matter perturbations as functions of real matter and induced 
metric perturbations 

5p = 5p(M) + ^(6-f^ + -^)-24^j , 

/«(4) I « J (b) 

5t, = <5vr(M) + 3 [ • 

Substituting the above expressions for the effective matter perturbations 5p, 5q, 5P and 
Sti in equations ( p^|l02|) , and replacing 5t:(^m) by zero, one obtains for a perfect fluid on a 
DGP brane 



^(''^ = 6^ { ^^^(0.0) + 0) + ^^?,o) + ^^^'C } ^ , (114) 



L " ) {h) 

«(4)'^P(M) = I -^'^ACCO) + ;^^'^A2(0,0) + -^^'^(1,0) + ^^^(2,0) + ;^^'^A(2,0) 
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+^^'C^o,i) + ^'C(i.i)}(,)' (116) 
k2 - i^^rP ^HQCP ^^^rP ^orP ^ rP 

+^^{3,0) + 'jp^iifi) + ^^'^(0,1) + ^'(^a,!) + ' (11'^) 

where the coefficients C are given in appendix and our normalization is chosen so that 
they are dimensionless. As for the RS model, we obtain from those equations a boundary 
condition for Pa expressible as 

FDGPmib) + GDGPmib) = 0, (118) 

where F£)qp and Gogp s-i'g polynomials of the cosmic time-derivative, 9^, with time 
dependent coefficients known from the background solution. 

The boundary conditions ( p.l2D and ( |118| ) have an unconventional form, since they 
involve derivatives of Vt and Vt' along the brane. However one can easily see, in the simple 
case where the master equation in the bulk is a wave equation, that such boundary con- 
ditions lead in this case to a well posed problem. This is shown explicitely in appendix 
We expect this to be true as well here, since those boundary conditions and Mukohyama's 
master equation encodes a priory all the physics of the problem. This can indeed be ver- 



ified numerically. Equations ([29| ) and (112) (respectively (118)) should be then all what 
is needed to solve for the evolution of RS (respectively DGP) brane world cosmological 
perturbations once initial conditions are supplied in the bulk as explained in subsection 



2.2 (see also appendix 0). Those two equations play for the brane world cosmological 



perturbations an equivalent role to the one played by equation (|106| ) for 4D adiabatic cos- 
mological perturbations of a perfect fluid. Once O is known in the vicinity of the brane, 
one can compute all the induced metric and real matter perturbations from equations 
([10^111 ) (respectively ( pl3HlT7| )). We stress again here that one can verify that the real 



matter, induced metric perturbations defined from those equations in terms of the master 
variable, as well as the Weyl's fluid perturbations given by equations (|95|), ( |96| ) and (97) 



verify together identically (i.e. whatever the function Q) the perturbed effective Einstein's 



equations (11) and conservation equations ( pi| ) and ( p^ (which components are gathered 
in appendix |B.2| ). As a last remark, we would like to comment on the GN coordinate 
system we used to derive the boundary condition. In general, there is no warranty that 
this coordinate system will not break down at some distance in the bulk. However, as we 
have seen above (from equation (^2|)) 17 is a gauge invariant quantity. Since the bound- 



ary condition ( |112D (or ( 118D ) is here expressed solely in terms of the master variable, it 



keeps the same form in a non GN coordinate system and can always be imposed on the 
background brane trajectory. 
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A Some results for the background metric 

In this appendix we give some results for the background bulk and brane metric. 
A.l Bulk relations 

The Einstein's equations (^) together with the cosmological ansatz ( |T8|) lead to the fol- 
lowing relations valid in the bulk (for a spatially flat metric on the brane) 
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The explicit form of the metric in the GN coordinate system ( [LSD can be obtained from 
those equations [37|. We wih not reproduce here the general results, but give here first 
the metric for the case of a Minkowskian bulk (with vanishing bulk cosmological constant 
and Weyl's tensor). In this case one finds (for a brane with positive effective tension) 



a 
n 



1 - \y\a{b)/a{b)- 



(123) 



In the above equation the unknown function of cosmic time aj-^) has to be determined 
solving the brane Priedmann's equations which are given in subsection A. 2. In the more 
general case where one only assumes a vanishing bulk Weyl's tensor, one gets the following 
relations valid everywhere in the bulk 
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A. 2 Brane relations 



Using the Junction conditions ( pj| ) and ( pOD a nd solving for the bulk metric one obtains 
the following brane Friedmann's equations |37| 
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p{2p + 3P), 



together with the conservation equation (^2|). The relations between effective matter 
energy density p and pressure P, entering into those equations, and real matter energy 
density and pressure P(m) depend on the model considered. In the case of the RS 
model (model(I)) one gets from equation (^) 
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P(M) + A(4), 
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For the DGP model (model (II)) one has ||T 



3 

P = PiM)--2--2^ (128) 



P = + 72+2- • (129) 



Using the above formulas, one can then solve for the explicit time dependence of the metric 
(see 13 for the RS model, and ll|,||] for the DGP model). 



B Scalar cosmological perturbations in 4D GR and on the 
brane 

In this appendix we first gather some results on standard 4D scalar cosmological pertur- 
bations and then give the analogous equations for brane worlds, which can be derived 



linearizing the effective Einstein's equations on the brane (11) and conservation equations 
and (0) (see e.g. H). 

B.l 4D cosmological perturbations 

We gather here some useful results concerning standard 4D scalar cosmological perturba- 



tions, we refer to the numerous review articles on the subject for more details (see e.g. [38|) 
We work in the longitudinal gauge, where the (4D) perturbed line element reads 



ds^ = -(1 + 2$)dt^ + a^(t)(l - 2^)6ijdx'dx^, (130) 

with a{t) being the background scale factor. The matter energy momentum tensor T^*^^ 
is given in the background by 

ri^^)'^ = diag(-p(M) , P(M) , P(M) , P(M) ) , 
while the perturbed matter energy momentum tensor is defined as in equations (41-43). 



A scalar gauge transformation given by 

t t + St, (131) 

x' x' + 5'^5xy, (132) 

induces the following transformation on the perturbed energy momentum tensor 

SP{M) ^P{M) - P{M)^t, 

^P(M) ^P(M) - P(M)5t-, 

^q(M) ^q(M) + {P{M) + P{M)) St, 

This means in particular that the perfect fluid condition ( |103| ) and the adiabatic condition 
( |104| ) are gauge invariant statements. The perturbed 4D Einstein's equations 

4)<^rif^=^GW, (133) 
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then lead to the four scalar equations 



2- „a/^a^\2 



4)5/.(M) = -6-[^ + -'^j+^A^, (134) 
^IMm) = -24>-2^^, (135) 

2 ^' — $ 

K(4)'5^'(A/) = 2 (2- + 4) ^ + 2-^ + 2^' + 6-^' + -^A($-^). (137) 

\ Qj Qj J Qj CL oQi 

Those equatfons are respectively the 00, Oi, traceless part of ij, and trace of ij components 
of the scalar part of Einstein's equations (|133| ). While the perturbed energy-momentum 
conservation equation 



leads to 



I)^T^f ) = 0, 



SP{M) + 3^(5/3(M) + SP{M)) + ^^Q{M) - ^{P{M) + P(A/))^ = 0, (138) 

2 

5q(^M) + ^^^QiM) + ^P{M) + (PiM) + P{M)) ^ + 3 A(5^(M) = 0. (139) 



Those two last equations can also be derived from the previous four ( ^3ijT37|) , thanks to 
Bianchi identities. 

B.2 Cosmological perturbations on the brane 

The perturbed 4D effective Einstein's equations (pTj) read 



with the scalar perturbations of 11 given by [20 



6U] = S]5Ut + ( a;. - -6}A ] 6Utf, 



1 



with 



OTt = l{{p + P)6p + p5P) 
OTtf = -^{p + 3P)6tt 



<5n° = lp6q\,. 



With the notations of the text, this leads to the following perturbed effective Einstein's 
equations on the brane 

= -4,^.,.,-{6^(* + f*)-lA*}^^^, (140) 
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^(5) 

6 



p6q = -4),5g(^)-|2vi/ + 2-cI>j , (141) 

L ) (b) 

'^U) o r ^> - ^ 1 

-^{3P + p)5n = K^4)57r(£) + |— ^1^^^, (142) 

'^i(p + P)5p + p5P) = -k2,)5P(^) + |(4- + 24U + 2-6 + 2^' 

6-^' + -^A($-^')l . (143) 
a 6a^ J (fc) 

Those equations are respectively the 00, Oi, traceless part of ij, and trace of ij com- 
ponents of the scalar part of the effective Einstein's equations (0). While the effective 
matter energy-momentum conservation equation ([l^) and Weyl matter energy-momentum 
conservation equation (0) lead to 



{5p + 3-{5p + 5P) + ^6q-3iP + p)4>\ = 0, (144) 

(5q + 3-5q + 6P+{P + p)<S> + -A5TT\ = 0, (145) 
la 3 J (5) 

|V(£) + 3^(<5p(^)+<5P(£))+^J<7(^)|^^^ = 0, (146) 



(3-5q + Ad7r-6p]] . (147) 
\ a J } (b) 



Again, those four equations are not independent from the four previous ones thanks to 
Bianchi identities. 

C Expression of the coefficients C 

Here we give the expressions for the coefficients appearing in equations ( 113| -117). 



Cl(o,o) = 1 + /t(T + 1), 

C(to) = 3{1 + /t(T + 1)}, 

Cdo) = -3/t(1 + T), 

Cfo.i) = 3{/t(T + 2)(T + 1)-T}, 

^A(0,0) = 1 + /x, 

C(to) = 3(1 + /t), 

C*i) = 3/x(2 + T), 
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-3(8 + 6T + r'^)H^HH + (4 + T)HH'^H + (2 + T)H'^{{A2 + 19T + T'^)H'^ 
+(2 + T)H-) + if2((54 + 30T + AT"^)!!^ - 2(2 + T)^^ + (2 + T)HH- ))\ , 
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eg,!) = 3/x(2 + T)2, 
where we have defined /t by 
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D Initial data and well-poseness 

We discuss here some issues related with the specification of initial data and the well- 
poseness of the differential problem defined in sections ^ and We will assume that 
initial data in the bulk are provided as indicated in the latter section and that a boundary 
condition on the brane is given by equation ( |112| ) or ( |118D of the former. We first discuss 
the relation between the bulk initial data and initial data on the brane for real matter and 
induced metric perturbations. The peculiarity of the boundary condition ( |112[ ) (or ( |118D ) 
is that this boundary condition involves derivatives of the master variable (and of its 
normal derivative 17') along the brane. In order to shed light on this type of boundary 
condition, we discuss in subsection p.2| a similar case where the PDE in the bulk is simply 
a wave equation. 

D. l Initial data 

Let us first consider the case where the initial data in the bulk are specified on a curve 
which is non characteristic in the vicinity of the initial event A chosen on the brane. 
This is the case for the initial curve (AD) or (AEF) on figure || (we refer here to figure 
1^ which concerns the DGP model, however the discussion also applies for the RS model). 
In this case we specify initial data in the vicinity of A by giving Cauchy type of data. 

E. g. on (AE) we specify 0, and its derivative normal to (AE). From these bulk data, 
one can get in A all the partial derivatives of Q of order greater or equal to one (this 
is because the initial data curve is non characteristic in A). This implies that the values 
of O, Q', and all their derivatives with respect to cosmic time are known in A from the 
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initial data in the bulk. These value should be compatible with the boundary condition 
(|112| ) (or ( |118D ). Furthermore this also implies that the values of the real matter and 
induced metric perturbations <5p(A4')5 '5-P(a/); '^'i'(M); ^(6) ^^'^ known initially in A from 
equations ([lO^|TTT|) or (|ll3| - p^ . 

Let us now turn to discuss the case where the initial curve is characteristic in the 
vicinity of A (The problem has then some similarities with the so-called Goursat problem) . 
This is for exemple the case if we specify initial data along the curve (AC) of figure |l|. Since 
this is a characteristic curve we choose only O along it, and can not specify freely there the 
normal derivative of fi. This implies that the derivatives of with respect to y and t of 
order greater or equal to one are not determined in A from bulk initial data. This is also 
true for the initial values of the real matter and induced metric perturbations. However, 
one can specify in A, in addition to characteristic initial data along (AC) (which allow to 
get the values of 9"0/ar", with n > 1, in A), the derivatives d'^Vt/dX" (with n > 1). If 
we do so, 0, Q' , as well as their derivatives with respect to cosmic time along the brane, 
will be known in A. This is also true for the initial values of real matter and induced metric 
perturbations. It is interesting to note that we have here the same freedom, by specifying 
X— derivatives in A, to specify initial data for real matter and induced metric as in the 
same case (i.e. adiabatic perturbations of a perfect fluid) in standard 4D cosmological 
perturbation theory. Namely the specification in A of and ^'(b) is enough to obtain 
the initial values of the other perturbations. 



D.2 Well-poseness for the wave equation in the bulk 

In the simple case where the master equation is the wave equation 

dxdrn = 0, (148) 

here written in its characterstics coordinates, it is simple to see that boundary conditions 
of the form ( |112| ) or ( [118| ) leads to a well posed problem, once initial data are chosen as 
explained above. We only discuss here the case where the curve on which we give initial 
data in the bulk is initially non characteristic. We will for simplicity assume that the initial 
data curve in the bulk is the curve (AEF) of figure |l|. Let us then consider a given point 
P (which is in the future of A) along the brane trajectory in the two dimensional space 
{X,T). We call Xp and Tp its coordinates. From the equation ( |148[) we know that 
is constant along the characteristic T = Tp (a similar discussion would hold for differential 
operators for which one can build Riemann invariant along T = constant characteristics). 
This means that (9t^ is known at P from the knowledge of the similar quantity at the 
intersection point Q of (AEF) and the T = Tp characteristic (the latter is given by initial 
data). From this we know then fi' in P as a function of Q, the brane trajectory and the 
initial data in Q. Inserting this expression of in the boundary condition ( [112| ) or ( p. 181 ) 
leads to an ordinary differential equation along the brane for J7. The initial time derivatives 
of $7 needed to solve this equation are known from bulk initial data, as is explained above. 
Then, solving for Q along the brane, we get a boundary condition along the brane which 
takes a familiar form. This enables to conclude that the problem is indeed well posed in 
the case considered here. 
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